Hard Discs on the Hyperbolic Plane 
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We examine a simple hard disc fluid with no long range interactions on the two dimensional space 
of constant negative Gaussian curvature, the hyperbolic plane. This geometry provides a natural 
mechanism by which global crystalline order is frustrated, allowing us to construct a tractable model 
of disordered monodisperse hard discs. We extend free area theory and the virial expansion to this 
regime, deriving the equation of state for the system, and compare its predictions with simulation 
near an isostatic packing in the curved space. 



The precise connections between the virial equation 
of state, the crystalline transition, glass transitions, and 
random close packing remain elusive P, 0] • The system- 
atic calculation of the equation of state is perturbative 
and analytic in volume fraction, </> [j| , and thus it is diffi- 
cult to comprehend how it can, taken at face value, probe 
the nonanalytic, large volume fraction regime where crys- 
tallization occurs. At the other extreme, crystal free en- 
ergies are well modeled by Kirkwood's free volume theory 
[3], especially near their maximum density. How can we 
probe, in three-dimensions for instance, the crystalline 
phase transition in hard spheres at <f>x ~ 0.5 and the 
possible existence of a maximally-random jammed state 
near 4>rcp ~ 0.64 [1, 5]? Focussing on two-dimensional 
systems as a simple starting point is fraught with diffi- 
culties. Unlike the situation in three dimensions where 
local tetrahedral and icosahedral packing competes with 
global fee crystallization it is problematic to prevent 
crystallization in simulations of monodisperse discs Q , as 
triangular close-packing is commensurate with crystalline 
order. Here we formulate the problem on a curved sur- 
face which serves to frustrate global crystalline order in a 
manner analogous to that in three dimensions. In partic- 
ular, we study hard discs on the hyperbolic plane, H 2 , at 
a curvature near a known regular tesselation. This model 
is numerically tractable and we find the equation of state 
for hard discs via molecular dynamics. We compare our 
results to free-area theory for the packing derived from 
the nearby tesselation. 

Experimentally, one could consider crystalline order or 
packing problems on any of the many bicontinuous struc- 
tures in lipid and diblock phaes Q and present in the cell, 
such as the periodic cubic membrane in the mitochondria 
of Chaos carolinensis ||. These bicontinous structures, 
often close to minimal surfaces, have negative intrinsic 
curvature and the regular tesselations we will consider 
here should serve as a scaffolding to model transitions on 
these surfaces. 

On a two-dimensional manifold with constant intrin- 
sic curvature K, the area of a circle of radius r is 
a = A(r) = -1-kK- 1 [cosh(y/^Kr) - l] (note that the 
K — > limit recapitulates the Euclidean plane). When 
we insert a disc of radius r onto our manifold, it excludes 
an area of radius 2r to the other discs. Thus the num- 



ber of discs that cannot be inserted once the first disc is 
placed A(2r)/A(r) ~ exp{ \J— Kr} grows exponentially 
with curvature and so we should expect that large, neg- 
ative curvatures will lead to highly correlated disc pack- 
ings due to the harsh entropic cost of this excluded area. 
Equivalently, the kissing number, nkiss, the number of 
spheres which may simultaneously touch a central sphere, 
as a function of the curvature, K, is [lo| : 
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To explore this systematically, we study the virial expan- 
sion for the pressure P in terms of the number density p, 
written P = kBTJ2^=iB n p n where the B\ term arises 
from the usual ideal gas law. While it may be tempt- 
ing to assume B\ = 1, the ideal gas on a curved man- 
ifold requires some thought; the circumference of a cir- 
cle, C(r) = 2ir sinh(\/— Kr) /yf—K, grows as fast as A(r) 
making the infinite area, thermodynamic limit problem- 
atic. 

The textbook approach to the ideal gas law begins with 
particle-in-a-box eigenenergies and takes the large area 
limit to calculate the partition function for a single par- 
ticle, Zi = V/X^h where A T = ^2nh 2 /(mT) is the ther- 
mal wavelength. Fortuitously, the spectrum of the Lapla- 
cian on a general Riemann surface has been well-studied. 
McKean and Singer ll| developed the Weyl expansion 
12] for the partition function of a single particle on a 
periodic two-dimensional domain D: 
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where the sign of the first correction depends on Neu- 
mann (+) or Dirichlet (-) boundary conditions and H is 
the mean curvature. Now we see that the circumference 
directly affects B\ and the only way to avoid the problem 
with the boundary is to consider periodic boundary con- 
ditions, i.e. no boundary. We note that for hard discs the 
virial expansion only includes local interactions among 
clusters of particles. As long as the clusters in question 
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do not wrap around the periodic space, we ignore the 
periodicity and perform the cluster integrals on H 2 . In 
order to ensure that this is possible, we may continue 
to increase the area available to a cluster by constructing 
ever higher genus manifolds [l5| and, thus, we can neglect 
"winding" of the discs around the periods. For manifolds 
with constant K < 0, Zi = A[X T 2 + K/(12w)] > when 
Ay < — 12ir/K. In other words, we are forced to work 
in the regime where the curvature of the manifold does 
not probe the quantum regime of the gas. In this regime, 
the pressure P = NttBTdlnZ/dA is independent of the 
proportionality constant. 

To calculate the virial coefficients, we employ the 
Poincare disc model of H 2 , popularized by Escher [lij ]. 
which conformally maps the hyperbolic plane with cur- 
vature K to the disc of radius \J — K~ Y on K 2 . Using the 
coordinates (x, y), the metric is [l5[ 
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ds 2 = 



4 (dx 2 + dy 2 ) 
[1 + K (x 2 + y 2 )f 



(3) 



where x 2 + y 2 < —1/K. Geodesies in this model are 
arcs of circles which intersect the bounding circle nor- 
mally; circles in the hyperbolic plane remain circles in 
this model, but their coordinate centers and radii vary. 
The radius of curvature sets a lengthscale and, conse- 
quently, any system on H 2 is manifestly not scale invari- 
ant. 

We take advantage of this lack of scale invariance to 
probe higher curvatures not by changing the ambient 
space, but simply by considering larger and larger discs. 
We consider discs of radius rv 1 —K. Without loss of gen- 
erality we fix the Gaussian curvature to be if = — 1. 
Dimensional analysis may always be employed to insert 
the appropriate factors of K . As in flat space, the sec- 
ond Virial coefficient, B%{r), is given by half the ex- 
cluded area of a single disc: B 2 {r) = 7r [cosh(2r) — 1]. 
As r — > 0, B 2 ~ 27rr 2 , we recover the d = 2 Euclidean 
result: B 2 (r) — (4-7r) 2 r d /r(|) d. Though it may been 
tempting to associate the higher connectivity at larger 
curvatures with higher dimensions, we know of no natu- 
ral identification. 

We turn to numerical integration to evaluate the higher 
virial coeffecients. In addition, as is traditional in flat 
space [13] , the higher coefficients B n are reported in units 
of a n_1 = [A(r)]" _ to suppress size effects and generate 
the expansion for Pj '(pfc^T) in the area fraction, <f> = pa: 
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We have checked the r — ► limit in all cases, and hap- 
pily find agreement with the known Euclidean results. 
We plot B 2 /a, B 3 /a 2 , -B4/a 3 , and B 5 /a 4 in Figure Q] as 
a function of r. To the accuracy of our calculation, B% 
and £>4 are positive, while B§ may have a zero crossing. 
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FIG. 1: (Color online). The second, third, fourth and fifth 
virial coefficients in appropriate units of the area, a, as func- 
tions of disc size (for K — —1) r. 



There is evidence that, in higher dimensional flat spaces, 
the virial coefficients alternate in sign 17, 18| and hence 
the leading singularity that controls the radius of con- 
vergence of the expansion sits on the negative real axis 
[l8| . The location of this singularity, and hence the range 
of applicability of the virial expansion itself, remains an 
open question in D = 2, 3. Whether or not B n (r) crosses 
zero for the negatively curved system at any value of r is 
the subject of future work. 

Since our confidence in the reliability of the virial ex- 
pansion wanes with ever larger curvature, we are forced 
to look at highly correlated structures on H 2 . To start 
we consider regular packings on the hyperbolic plane. We 
begin with a regular lattice of Voronoi cells and consider 
a tesselation of F p-gons, with q meeting at each vertex. 
The total number of vertices is pF/q, the total number 
of edges is pF/2, and the total number of faces is F. The 
Euler character \ of the surface constrains these three 
numbers and it follows that 



X = V- E + F=pF 
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For the flat periodic plane (i.e. the torus) x = 0, so 
p^ 1 + q^ 1 = 1/2 and we find the three tilings rep- 
resented by the Schla.fi symbol {p, q} — {6,3}, {3,6}, 
and {4,4}, hexagons, triangles, and squares, respec- 
tively. The spherical topology is the only one with pos- 
itive x = 2, admitting the five Platonic solids {p, q] = 
{3, 3}, {3, 4}, {3, 5}, {4, 3}, and {5,3}. Finally, we turn 
to the pertinent geometries with negative Gaussian cur- 



vature, for which x < and p 



< 



In this 



case there are an infinite number of integral pairs {p, q} 
and an infinite number of regular tesselations. Triangular 
packing around each vertex, q = 3, corresponds to close- 
packing. We focus instead on isostatic packings, those 
for which each disc has p — 4 nearest neighbors. Recall 
that an isostatic packing is one for which the number of 
force balance equations is equal to the number of degrees 
of freedom. For N frictionless particles in d dimensions 
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with coordination z, there are dN degrees of freedom and 
Nz/2 forces to balance: isostaticity requires z = 2d 

We find the area fraction of {p, q} through appli- 
cation of the Gauss Bonnet Theorem and hyperbolic 
trigonometry. When there are F identical polygons, 
the area of each can be found from the Euler charac- 
ter x- FA p - gon — J dA — — J KdA = —2irx, where we 
have used K = — 1. The radius of the incircle can be 
determined via the dual law of hyperbolic cosines [ioj ]. 
cosh(r) = cos(7r/g)/ sin(7r/f>). We find the area fraction 
for the {p, q} tesselation: 
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We recover the appropriate packing fractions, 7r/3v3, 
7r/4, and 7r/2v3, for the Euclidean tilings, {3,6}, {4,4}, 
and {6,3}, respectively. As p — > oo for close-packed 
q = 3 tesselations, we find max = 0.9549, the known 
packing fraction of the best packing on the hyperbolic 
plane at any curvature 2l|. The lowest packing fraction 
for an isostatic configuration (p = 4) is </>j SO = 0{4 iOO } = 
\/2-l» 0.4142, an area fraction far below that for the 
Euclidean square lattice 7r/4 w 0.785. Also of note is the 
broadening of the range in <p that supports stable con- 
figurations: </>max — 4>\so ~ 0.54 on H 2 , as compared to 
a range of roughly 0.13 in flat space. This is consistent 
with our finding that the virial expansion breaks down 
at lower area fractions on H 2 . Indeed, the {4, oo} tes- 
selation requires cosh(r*) = \f2 so that B2(r*)/A(r*) = 
i(2)/(v2— 1) = l/0{4,oo}- Though possibly a mathemat- 
ical profundity, the fact that i?2(r*)(/>{4 i00 }/A(r*) = B\ 
means that the first correction to the ideal gas is precisely 
equal to the non-interacting result and strengthens the 
case that we must study the large <f> regime. 

The presence of these large corrections, even at low 
area fraction, necessitates an approach to the high- 
density regime. We turn to free-area theory, known to 
be exact near close-packing [22,[23J. In flat space, scale 
invariance ensures that the shape of the free-area region, 
available to each disc's center of mass, has precisely the 
same shape as the corresponding Voronoi cell. However, 
the rigidity of H 2 forces us to consider how the shape of 
the cell changes as the size of the discs change. In Figure 
[5] we show how we break up the Voronoi cell to calculate 
the resulting free area. 

The free-area can be found by starting with A p _ gon 
and subtracting the area of each of the p strips, ^4 s trip = 
(s — 28) sinh(r) [24| and each of the p corners A comer = 
7r — a — 2n/q. The hyperbolic law of cosines [2(| allows us 
to find cosh 8 — cos(a/2) csc(7r/g), coshs = [cos(27r/p) + 
cos 2 (7r/g)] csc 2 (7r/g), and sin(a/2) = cos(ir/q)/ coshr. 
We can relate r to </>, p, and q and write j4f r0 c in terms 
of p, q, and 4>: setting r = 1 - *2 and r pq = 1 - x4>i p q} , 



FIG. 2: Partition of a {p, q} Voronoi cell into its free area cell, 
p corners, and p strips. Note that the boundaries of the free 
area cell are not geodesies. 
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Unfortunately, this expression is cumbersome to the 
point of uselessness, so we instead turn to Figure^ where 
we plot the free-area equations of state for various {p, q} 
tesselations using Pfa = —Ta^ 1 (j) 2 -^ln(Af rcc ). Note 
that a is set by A p _ gon (j) = a. 

Finally, we turn to disordered packings. We employed 
constant energy molecular dynamics on a periodic hy- 
perbolic octagon with sides identified to give the topol- 
ogy of the 2-torus [15j. Note that we cannot vary by 



changing the disc size, as we must in free area theory, 
since that is, by the absence of scale invariance, equiva- 
lent to a change in the background curvature. Instead, 
we change <fi by changing the number of discs in the pe- 
riodic octagon. We could choose a larger, higher genus 
periodic cell on which to do our simulation in which case 
we could probe a larger number of area fractions, filling 
in the data points in Figure [31 though this is compu- 
tationally intensive. Although the simulation and free- 
area theory share a common endpoint, the pressure in 
the latter samples ever-changing curvatures on its way 
to close packing. We chose to study curvatures r = 0.53 
and r — 0.531, just below and just above that for the 
{4,5} tesselation, r{4 j5 } = 0.5306. At these curvatures 
no regular tesselation is allowed, and so non-crystalline 
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FIG. 3: (Color online). The equation of state for free area 
theory for a range of {p, q} tesselations. Note that the pres- 
sure curves diverge at their maximum packing fractions and 
that they differ even at low volume fractions, (f>. The verti- 
cal dashed lines indicate the maximum volume fraction for 
the tesselations listed on top of the figure. The inset com- 
pares our numerical results with the {4, 5} free-area theory at 
larger values of the pressure. The larger, light blue circles are 
the result of our molecular dynamics simulation for curvature 
just below the {4, 5} tesselation and the smaller, black circles 
are for curvature just above {4,5}. The number of particles 
ranges from 1 to 9. 



arrangements give the best packings [25| . Thus we expect 
to probe disordered configurations even at high densities 
with this choice. Indeed, the simulation suggests that the 
disordered system attains higher pressures less rapidly 
than the nearby, ordered, isostatic configurations, par- 
ticularly at larger values of </5 (Figure [3]) . The r — > 
limit should be a way to study these disordered packings 
on the Euclidean plane. Fortuitously, roughly 10 discs 
fill the periodic octagon, the smallest periodic region on 
H 2 , and thus we could quickly approach a dense pack- 
ing. We require many more discs to study the Euclidean 
limit; this work is in progress. 

We studied a hard disc gas through simulation, at cur- 
vatures close to that for a regular tesselation and com- 
pared with free-area theory. Thus we have a model which 
approaches disordered packings of monodisperse discs for 
which we can understand the equation of state. Note that 
our approach allows us to ponder the structure of the 
virial expansion: in d dimensions, the 2d th virial coeffi- 
cient, B2d is the required to probe isostaticity, while crys- 
tallization requires the n^ ss virial coefficient. Though in- 
creasing the dimension d pushes nkiss ever larger [26j , the 
study of disordered packings forces us to calculate ever 
more virial coefficients. The hyperbolic plane sidesteps 
this problem by keeping isostaticity at p — 4. Further 
work will focus on other tesselations, in particular, the 
Euclidean limit, and on mapping the results and predic- 



tions of H 2 onto the bicontinuous structures present in 
biological systems. 
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